We claim: 

1 1 . A computer-assisted method for creating and displaying a model of a molecule in which 

2 the residues of said molecule that are affected by binding of a ligand to said molecule are 

3 highlighted using a programmed computer including a processor , an input device, and an 

4 output device, comprising, 

5 inputting into said programmed computer, through the input device, data including an 

6 identity and three-dimensional coordinates of said molecule, 

7 generating, using the processor, a conformational ensemble for said molecule, 

8 wherein said conformational ensemble is comprised of a number of conformational states of 

9 said molecule sufficient to achieve convergence to a stable ensemble average, 

1 0 determining, using the processor, a predicted Gibbs energy for each of said 

1 1 conformational states, 

1:2 identifying, using the processor, binding competent conformational states, 

13 modifying said predicted Gibbs energy for said binding competent conformational 

14 states to form a modified Gibbs energy using an equation 

£5 

16 (l+Ka,i[X]) 

17 AGi = AGi° - RT In - 



(l+Ka, 0 [X]) 



2- 0 wherein AGi is said modified Gibbs energy, AGi 0 is said predicted Gibbs energy, R 

21 is a gas constant, T is an absolute temperature, Ka, 0 is a binding constant to a reference state, 

2 2 Ka„ is a binding constant to a given binding competent conformational state, and [X] is a 

2 3 concentration of said ligand, 

24 forming, using said processor, a first and second set of Gibbs energy values, wherein 

2 5 said first set of Gibbs energy values represents Gibbs energy values in the presence of said 

2 6 ligand and said second set of Gibbs energy values represents Gibbs energy values in the 

2 7 absence of said ligand, 

2 8 computing, using the processor, a first set of probabilities from said first set of Gibbs 

2 9 energy values and a second set second of probabilities from said second set of Gibbs energy 

3 0 values, wherein said first set of probabilities represents probability values in the presence of 
3 1 said ligand and said second set of probabilities represents probability values in the absence 
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of said ligand, 

calculating, using the processor, a first set of residue stability constants in the 
presence of ligand from said first set of probabilities and a second set of residue stability 
constants in the absence of ligand from said second set of probabilities, 

comparing, using the processor, said first and second sets of residue stability 
constants, 

selecting affected residues of said molecule, wherein said affected residues exhibit a 
difference in stability constant in said first and said second set of stability constants, 

creating, using the processor, a model of said molecule wherein said affected residues 
are highlighted, and 

displaying, using the output device, said model of said molecule wherein said 
affected residues are highlighted. 

Claim 2. The method of claim 1 wherein said step of identifying binding competent 
conformational states comprises the steps of 

ascertaining a contribution of each residue of said conformational states in said 
ensemble to the Gibbs energy of binding for said ligand, 

defining binding determinant residues, wherein said binding determinant residues 
contribute in excess of about 100 cal/mole to the Gibbs energy of binding, and 

establishing as binding competent conformational states those conformational states 
in which said binding determinant residues are in a native conformation. 

Claim 3. The method of claim 1 wherein said molecule is a protein. 

Claim 4. A computer program product comprising a computer-readable medium having 

stored thereon a computer program including instructions for causing a computer to: 

receive data including an identity and three-dimensional coordinates of a molecule, 
generate a conformational ensemble for said molecule, wherein said conformational 

ensemble is comprised of a plurality of conformational states of said molecule, 

determine a predicted Gibbs energy for each of said plurality of states in said 

conformational ensemble, 
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8 identify binding competent conformational states, 

9 modify said predicted Gibbs energy for said binding competent conformational states 
10 to form a modified Gibbs energy using an equation 



12 a * „ (l+Ka,i[X]) 

13 AGi = AGi - RT In - 



(l+Ka, 0 [X]) 



15 

16 wherein AGi is said modified Gibbs energy, AGi 0 is said predicted Gibbs energy, R 

17 is a gas constant, T is an absolute temperature, Ka, 0 is a binding constant to a reference state, 

1 8 Ka„ is a binding constant to a given binding competent conformational state, and [X] is a 

1 9 concentration of said ligand, 

2 0 compute a first set of probabilities from said first set of Gibbs energy values and a 

2-1 second set second of probabilities from said second set of Gibbs energy values, wherein said 

22 first set of probabilities represents probability values in the presence of said ligand and said 

2 3 second set of probabilities represents probability values in the absence of said ligand, 

24 calculate a first set of residue stability constants in the presence of ligand from said 

25 first set of probabilities and a second set of residue stability constants in the absence of 
2€ ligand from said second set of probabilities, 

2 ? compare said first and second sets of residue stability constants, 

28 select affected residues of said molecule, wherein said affected residues exhibit a 

2 9 difference in stability constant in said first and said second set of stability constants, 

3 0 create a model of said molecule wherein said affected residues are highlighted., and 
3 1 display said model of said molecule. 
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